The surface brightness of dark matter: 
unique signatures of neutralino annihilation in the Galactic halo. 



Carlos Calcaneo-Roldan* & Ben Moored 
Department of Physics, Durham University, Science laboratories, Durham, DH1 3LE, UK. 

(February 1, 2008) 

We use high resolution numerical simulations of the formation of cold dark matter halos to simulate 
the background of decay products from neutralino annihilation, such as gamma-rays or neutrinos. 
Halos are non-spherical, have steep singular density profiles and contain many thousands of sur- 
viving dark matter substructure clumps. This leads to several unique signatures in the gamma-ray 
background that may be confirmed or rejected by the next generation of gamma-ray experiments. 
Most importantly, the diffuse background is enhanced by over two orders of magnitude due to an- 
nihilation within substructure halos. The largest dark substructures are easily visibly above the 
background and may account for the unidentified EGRET sources. A deep strip survey of the 
gamma-ray background would allow the shape of the Galactic halo to be quantified. 
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I. INTRODUCTION 

Determining the nature of dark matter is of fundamen- 
tal importance to both Astronomy and Particle Physics. 
Both theory and observational data currently favour a 
universe with a matter density that is dominated by non- 
baryonic particles. Many candidates have been proposed: 
some are known to exist, others are more speculative 
(e.g., Ref. jl] and references therein). Structure forma- 
tion in a universe dominated by cold dark matter (CDM) 
has been extensively tested against observations and the 
model has proven highly successful at reproducing the 
large scale properties and distributions of galaxies ||[| . 
On the non-linear scales of galactic halos it remains to 
be confirmed wether the model can successfully repro- 
duce the observational data ||-|| . 

Direct detection in the laboratory is the ultimate tech- 
nique for verifying the existence of dark matter particles 
(see, Ref. §). However, even the most popular candi- 
date for dark matter, the neutralino, has a cross section 
that spans many orders of magnitude and the current lab- 
oratory searches are only just becoming sensitive to the 
cosmologically interesting parameter range. Presently, 
Astronomical observations provide the best insights into 
the nature of the dark matter, furthermore direct detec- 
tion relies on the existence of a smooth component of 
dark matter. 

Within the next few years indirect detection of neu- 
tralinos will provide interesting constraints on their possi- 
ble cross-section and masses. Neutralino-neutralino anni- 
hilation produces observable photons (as well as a host of 
other particles) that may be observed as a diffuse gamma- 
ray background from the halo surrounding the Milky Way 
as discussed in Refs. and more recently, in Refs. 
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Renewed interest in these predictions has recently 
arisen because of an unexplained component of diffuse 
high energy photons in the Egret data (e.g. Ref. p(i|), 
and also the possibility of an excess from the center of 
the Galaxy itself pi] , unexpected clumpy emission and 
the unresolved "discrete sources" |22j. Progress in this 
area will result from several new and sensitive gamma-ray 
surveys such as GLAST [|| and VERITAS @. 

The efficiency of the annihilation process is strongly 
dependent on the local density and the cross-section of 
the neutralino. Many authors have calculated the ex- 
pected flux from the Galactic halo using simple models 
for the expected mass distribution of neutralinos within 
the Galaxy (see, e.g., Refs. |l^,|l^|) or from its satellites 

EMU- 

Advances in computational Cosmology have lead to 
several recent breakthroughs that have direct relevance 
to the detection of dark matter. In particular, the numer- 
ical resolution that can be achieved using parallel com- 
putational techniques is now sufficient to study the in- 
ternal structure of dark matter halos that form within 
a cosmological context. The results of these simulations 
have important implications for indirect (and direct) de- 
tection of dark matter candidates. Most significantly for 
particle-particle annihilation, we are now confident that 
the central density profile of CDM halos follows a singu- 
lar power law down to small scales p5|-p8| . Thus we may 
expect a point like source of mono-chromatic gamma-rays 
emanating from the center of the Milky Way, where the 
annihilation rate will be very high. 

A second fundamental prediction of the CDM model is 
that previous generations of the merging hierarchy sur- 
vive within halos ps|] . Halos that accrete into larger 
systems may be tidally stripped of most of their mass, 
however their dense central cores survive and continue 
to orbit within the parent halos. This may present some 
problems for the CDM model since the predicted number 
of satellites within the Milky Way's halo is 50-100 times 
as many as observed pjj. If the CDM model is correct, 
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then only a fraction of these satellites must have formed 
stars and most of the substructure remains as dark ob- 
jects within the Galactic halo. 

The possibility of an enhanced gamma-ray back- 
ground from dark matter substructure was explored by 
Bergstrom et al. [P9| , who made simple assumptions as 
to the mean density and abundance of such clumps. We 
can now use the high resolution N-body simulations to 
directly measure these quantities. The simulations also 
allow us to study the influence of the halo shapes on the 
diffuse gamma-ray background and the intensity of the 
central halo emission that arises from the singular dark 
matter density profiles. This paper is organized as fol- 
lows. In Section || we explore the gamma-ray background 
that results from the smoothly distributed component of 
dark matter using both analytical and simulated halos. 
In Section HI we focus on the substructure within halos. 



Our conclusions are summarised in Section [V 



II. 



THE SKY DISTRIBUTION OF THE 
GAMMA-RAY BACKGROUND 



In what follows we will consider a flux of photons 
(or other particles) that are a by-product of the an- 
nihilation of dark matter particles within the smooth 
component of dark matter that surrounds the Galaxy. 
It is not our intention to discuss the details of neu- 
tralino interactions, a complete overview on these pro- 
cesses (and super-symmetric matter in general) can be 
found in Refs. |3Q,B1|. 



A. Model neutralino halos 

We calculate the gamma-ray flux along a given line of 
sight through a spherically symmetric galactic halo using: 
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f sight 
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where ip is the angle between the direction of galactic 
center and observation; p, the density of dark matter at 
distance I from the observer. We have summed up the 
dependence of the flux on neutralino mass and interaction 
cross section in the constant K . This is enough scope 
for the present discussion - it is straightforward to take 
our results and input a neutralino cross-section, < av >, 
and mass, M x to determine the absolute gamma-ray flux 
(where K is defined to be < av > /M 2 ). Our results can 
also be used to infer the sky distribution of other products 
of the annihilation, such as neutrinos or positrons. 

The line of sight distance, I, is related to the radial 
distance from the halo center, r, via 

r 2 = I 2 + R 2 a - 2lR a cos(V>) 



where R is our galacto-centric distance, taken here to 
have the IAU standard value of R Q = 8.5kpc and i[> 
is related to galactic coordinates (£, b) through 

cos(V') = cos(£) cos(b). 

For the halo density profile, p(r), we take the latest re- 
sults from the highest resolution numerical simulations of 
galactic halos carried out to-date || . These authors sim- 
ulated 6 different galactic mass halos with force resolu- 
tion of 0.5 kpc and mass resolution of 1O 6 M0. (Through- 
out the paper we will use the Hubble constant value of 
H a = 100/ikms~ 1 Mpc _1 and h = 0.5; as adopted for 
the simulations.) The best fitting density profile to these 
data is (subscript moore): 
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Where r is the distance from the halo center, a — 
f 200 1 Cmoore the scale radius for halos of mass « 1 X 
10 12 M Q . The virial radius of our fiducial Galactic halo, 
r 2 oo ~ 300 kpc, is defined as the radius of a sphere at 
which the mean overdensity is 200 times the cosmologi- 
cal mean density. (A central density profile of slope -1.5 
on galactic scales was also found by Jing & Suto |33j ] and 
confirmed as an asymptotic slope by Ghigna et al. ft-if.) 

We also compare this profile with that determined by 
Navarro et al. [p6| using a sequence of lower resolution 
studies (subscript nfw) (the main difference being that 
the central dark matter density profile has a slope of -1): 



Pnfw (?) 



Pnfv 



{r / a){\ + r / af 



(3) 



and the modified isothermal profile with a constant den- 
sity core (subscript is): 
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[1 + (r/a) 2 ] 3 / 2 



(4) 



The scale radius, a, is determined directly from the 
numerical simulations, except for the modified isother- 
mal model which we normalise to match the observa- 
tional rotation curve data (as in Ref. ]3q| ); ai S — 24.3 
kpc, a n f w = 27.7 kpc and a moore = 33.2 kpc (this 
radius is directly related to the concentration param- 
eter, c = r2oo/a). We normalise each density profile 
such that the peak circular velocity, v pea k = 200 km/s 
(the maximum of the v c = y/GM/r curve), which gives: 
p' is = 4.96 x 10 6 M Q kpc~ 3 , p' nfw = 5.11 x 1O 6 M kpc" 3 



and p' moore = 1.64 x 10 6 M Q kpc" 3 . We plot the effec- 
tive circular velocity profiles and density profiles of these 
model halos in Fig. |l|(a) and Fig. |l](b) respectively. 

In Fig. H we plot the flux, <p, along the line of sight 
through a spherical Milky Way halo using the above den- 
sity profiles as the observer looks towards the Galac- 
tic center at ip = 0°, to the Galactic anticenter at 
ijj = 180°. As expected, the central annihilation flux 
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depends strongly on the form of the inner density pro- 
file. At an angle of five degrees from the Galactic cen- 
ter, the ratio of fluxes from the three different profiles, 
moore:nfw:is are 1000:100:1. 

The peak central value depends upon the distance from 
the Galactic center that we are willing to consider inte- 
grating from - the flux slowly diverges for the density 
profile in Eq. (2). However, within a given radius, most of 
the neutralinos would have self annihilated leaving a tiny 
constant density core. We can estimate the size of this 
core using (nav)^ 1 = t h , where t h ~ 10 Gyrs is the Hub- 
ble time. Taking a typical cross section, av = 10~ 30 cm 3 
s _1 , and adopting the Moore et al. density profile we 
find that the annihilation radius within the Milky Way 
is approximately 4 x 10 - ' parsecs « 10 _12 r2oo- 

The total flux that arises within 5 degrees of the Galac- 
tic center using the Moore et al. density profile is a fac- 
tor of 20 larger than that found using the NFW profile 
(both integrated down to the annihilation radius calcu- 
lated above). 

B. Comparison with high resolution CDM 
simulations 

We can use the numerical simulations to compare di- 
rectly with the above predictions that were obtained 
assuming spherical symmetry. We refer the reader to 
Moore et al. (Ref. ||) for details of the numerical 
simulations .0 To construct the expected gamma-ray sky 
maps we choose a simulated dark matter halo at a red- 
shift z=0 that has a peak circular velocity of ~ 200kms _1 
and a total mass, within the virial radius, r2oo = 300kpc, 
of 1 x 10 12 M Q . This simulated halo is from the "Local 
Group" simulation and is close to our fiducial Milky Way 
cold dark matter halo that we adopted in the previous 
section. 

N-body simulations attempt to simulate a collisionless 
fluid of dark matter using discrete massive particles. We 
calculate the local density at the position of each particle 
by averaging over its nearest 64 neighbours. The observer 
is placed 8.5 kpc from the halo center (defined using the 
most bound particle in the simulation) and we sum up 
the flux of annihilation products along each line of sight 
using the discrete equivalent to Eq. ( ||): 

*& & ) = TT £ Pl^b)An(£,b) (5) 

LOS 

where £, b are galactic longitude and latitude respectively. 
The flux is binned in angular windows of size = 1° x 1° 
and in the radial direction in fixed increments Art = 
lkpc. 



The simulated dark matter halos are typically flattened 
oblate or prolate systems jjfj). We do not know a-priori 
in which axis the stellar disk would be located, therefore 
we show two all-sky maps using the same dark matter 
halo but viewed using two different locations for the "ob- 
server": Figure ||(a) and Fig. ||(b) has the observer lo- 
cated on the short and long axes respectively. Both of 
these plots show the enhanced brightening towards the 
halo center, as well as some clumpy substructure in the 
halo itself. Note that both the central halo and the cen- 
ters of the substructure halos are artificially "dimmed" 
in these plots due to the numerical resolution ~ 0.5 kpc, 
which sets a maximum density that can be resolved. The 
non-spherical shape of the halo is also clearly evident by 
inspecting the plots with different observer positions. 

Recent estimates for the shape of the Milky Way's halo 
(see, e.g., |$7) and references therein), suggest that it 
may be flattened with a short/long axis ratio of 0.5. An 
independent estimate from the orbit of the Sagittarius 
debris stars yields a nearly spherical dark matter halo 
p8| . The simulated halo that we have chosen to analyse 
represents a typical prolate CDM halo with a short to 
long axis ratio of 0.4, and intermediate to long axis ratio 
of 0.5. 

It is straightforward to estimate the effects of flattened 
dark matter halos by modifying Eq. ( 0) to accommodate 
triaxial shaped bodies. The simplest way to achieve this 
is to change from spherical coordinate r to 

<r2 _ x 2 +y 2 , 

? — b 2 ~r 7? 

where b > c for the oblate case, and b < c for prolate and 
we leave z as the axis of symmetry. A 2d visualization of 
these 3d shapes is illustrated in Fig. ||. 

In Fig. |5| we plot spherical, oblate (2:1) and prolate 
(2:1) versions of the integral in Eq. ( using the Moore 
et al. (1999) density profile. The observer is located on a 
plane parallel to the axis of symmetry, again at a distance 
R = 8.5kpc from the center of the halo. The halo shape 
leads to little difference towards the Galactic center, but 
at the anti-center prolate halos can be 100 times brighter 
than oblate halos. 

We can also compare the predicted angular flux with 
that measured directly from the N-body simulation. The 
annihilation flux is averaged in ten degree bins from the 
simulated dark matter halo, along a great circle from 
from the galactic center to its anti-center. This direct 
measurement of the flux is also plotted (as points) in 
Fig. g. These data are particularly noisy due to the large 
numbers of substructure clumps in the simulation - the 
spike at tp = 125° is due to a massive dark clump that 
happens to lie exactly along this chosen line of sight. 

III. SUBSTRUCTURE 



'Images, data and mov ies of these dark matter simulations 
can be downloaded from jittp: / /www. nbody.net. 
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A. Enhancement of global flux due to substructure 

Cold dark matter substructure clumps have singular 
density profiles that will be a significant source of annihi- 
lation products. The velocities and spatial distribution of 
dark matter substructure is unbiased with respect to the 
smooth dark matter background [Q. Therefore, to first 
order, substructure increases the global sky brightness 
in any given direction. However, the details depend on 
how much substructure survives within the solar radius 
and also on how far down the mass function substructure 
halos form and survive. 

First we will estimate the annihilation flux from clumps 
of dark matter that are known to exist in the Galactic 
halo i.e. the dark matter halos that surround the Magel- 
lanic Clouds and dwarf spheroidal galaxies. In fact, high- 
Energy gamma-ray emission from the Large Magellanic 
Cloud (LMC) was detected with EGRET by Sreekumar 
et al. |39[ in 1992 (although the origin of this emission 
was reported to be the interaction of cosmic rays with 
interstellar matter). 

We estimate the average flux, $av, from the dark 
matter halos that surround some of the principal struc- 
tures in the Local Group: The "Andromeda Galaxy"; 
M31 (vp ea k — 200kms~ 1 at a distance of 700 kpc), The 
Large and Small Magellanic Clouds (v pen k = 70kms _1 



and v, 



peak 



4Qkms -1 , respectively; both at a distance 
of 50 kpc), Draco (v pea k = lOkms -1 at a distance of 50 
kpc) and a small dark matter clump (v pea k — 2kms~ 1 at 
a distance of 10 kpc). A sketch of the geometry is given 
in Fig. | 

The total flux from a substructure halo at distance R c 
from the observer is 



®tot{Rc) = ^J P 2 (r)r 



'dr. 



(0) 



by considering the central A£l = 1° x 1° patch over each 
clump, we define the maximum integration limit in |6| 
above and the average flux is then 
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(we set AVI in steradians, so we may compare directly 
with the smooth flux of Section ||) . 

For the dark matter distribution within the substruc- 
ture clumps we use the Moore et al. profile, which pro- 
vides a good fit to the smallest, well-resolved substruc- 
ture halos. The concentration of CDM halos is a function 
of mass [Eo| and for the density profile in Eq. || this can 
be written: 



0.084 



(8) 



This defines the scale radius of each substructure 
clump: a M 3i = 33.3kpc, a LM c = 6.7kpc, a S MC = 

3.1kpC, dDraco = 0.5kpC, axiny = 0.05kpc. 



The integral in Eq. || diverges as r — > for the den- 
sity profile that we are using, however, even the smallest 
substructure halos will have a maximum density set by 
the radius within which most of the neutralinos would 
have self annihilated. We therefore present results for 
the average flux from these clumps as a function of the 
minimum integration radius R m i n /a in Fig. |7[ where a is 
the scale radius as defined above. 

For comparison, we plot the range of background emis- 
sion at the Galactic anti-center as the shaded line in 
Fig. [t]. The Tiny clump is only marginally visible above 
the background flux (depending on whether or not the 
Galactic halo is prolate or oblate) whereas most of the 
subhalos are easily visible. Also for comparison we have 
plotted the flux from the inner region of the galaxy which 
is the brightest of these sources. 

Although the Galactic halo is expected to contain just 
a few clumps more massive that the Magellanic Clouds, 
there are many thousands of smaller mass objects. The 
mass function of substructure is a power law close to 
dn(m)/dm oc m~ 19 or in terms of circular velocity 
dn(v c )/dv c oc v~ 3 ' 8 ( fl34|). Above a circular velocity 
Vpeak = 10 km s _1 and 1 km s _1 we expect the galactic 
halo to host roughly 1000 and 5 x 10 5 substructure halos 
respectively. Future simulations should be able to mea- 
sure how far down the mass function substructure halos 
can survive as well as to determine their central density 
profiles. (We note that the highest resolution simulation 
to date resolved the substructure within a dark matter 
mini- halo of mass 10 7 M Q . The force resolution was 10 
parsecs and the mass resolution was 10M Q allowing sub- 
structure with peak circular velocities as low as a few 
hundred meters per second to be resolved. The survival 
of substructure continues even down to this scale, where 
the slope of the power spectrum is close to -3.) 

We calculate the total flux from substructure using 
Monte-Carlo techniques. First we generate a list of peak 
circular velocities and positions of 5 x 10 5 substructure 
halos in the range of l-70km s" 1 . (Distances are ran- 
domly selected using the Moore et al. density profile 
and peak circular velocities are randomly assigned from 
a power law distribution scaling as t>~ 3 ' 8 ). For each lump, 
we estimate its total flux as in the previous cases, inte- 
grating Eq. ^| with a density profile scaled according to 
Eq. ^ for the concentration. 

In the absence of further constraints on the possi- 
ble value for R m i n , we use the same criteria as before 
and choose it to be a fixed fraction of the virial radius, 
Rmin ~ 10~ 12 r 2 oo- This corresponds to a mean density 
of w 10 22 M Q kpc -3 for the galactic halo. The results are 
not too sensitive to the value of the minimum integration 
radius as is apparent from inspection of Fig. (7| The total 
flux is then averaged over the entire sky and we repeat 
this process in order to estimate the variance. The cu- 
mulative distribution of flux above a given substructure 
peak circular velocity (S$tot) is plotted for ten of these 
random halo realisations in Fig. ^. 

It is evident from this plot that the effects of includ- 
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ing the entire mass spectrum of substructure is quite 
dramatic and boosts the expected flux from the smooth 
halo by several orders of magnitude. However, most of 
the flux arises from the subhalos with circular velocities 
larger than 10 km s _1 . Extrapolating to very small halos 
would not change the total flux by a large factor. 

To quantify the brightening of the background due to 
substructure, we have calculate the average flux due to 
all clumps with v pea k > 1 km s _1 within a spherical halo. 
The point plotted in Fig. [7] represents this contribution 
to the flux, where the error-bar is the la variation among 
the different Monte-Carlo models. From this plot we see 
that the flux due to substructure is over two orders of 
magnitude brighter than the smooth background from 
a spherical halo. We note that one needs to observe a 
fairly large fraction of the sky (> 100 square degrees) to 
ensure a significant number of clumps lie in the field of 
view. (Also note that the variance at high peak circu- 
lar velocities is due to the proximity of the largest few 
dark matter substructures, however, the mean total flux 
converges to similar values for each Monte-Carlo model.) 

B. The flux due to substructure in prolate and 
oblate halos 

Not only is the mean flux at a given position on the 
sky dominated by substructure halos, the spatial distri- 
bution of flux across the sky will be determined by the 
substructure. The convergence study by Ghigna et al. 
(Ref. Q ) shows that substructure halos trace the global 
mass distribution of the halo. Therefore, we can use the 
N-body simulations to generate Monte-Carlo distribu- 
tions of substructure halos and construct all-sky maps of 
the expected gamma-ray flux. We take a random particle 
from the simulation and assign a circular velocity from 
a distribution dn(v c )/dv c oc t>~ 3 ' 8 . For each sub-halo we 
calculate its total annihilation flux and then repeat the 
process until we have 500,000 halos above a circular ve- 
locity of lkms -1 . 

Figure ^|(a) and Fig. [|(b) show the resulting sky distri- 
bution of flux from sub-halos binned in one degree bins 
where the observer has been placed in the short and long 
axis of the simulation respectively. Large substructure 
halos, such as the Magellanic Clouds in our own halo, will 
contain its own gravitationally bound sub-halos which 
leads to clustering of gamma-ray emission in the all sky 
maps. 

Future observations may only be able to make deep 
strip maps therefore in Fig. [lo] we have binned the flux 
along lines of constant galactic I and b with the observer 
placed in the short and long axis of the global density dis- 
tribution. From these plots we can see that the emission 
from substructure peaks at the galactic center, as one 
would expect, this effect is not that different for spheri- 
cal halos than for prolate or oblate halos. 

Since the substructure traces the global mass distribu- 
tion, a prolate halo would also have a prolate distribution 



of satellites. Therefore we can study the variation of flux 
within smooth prolate or oblate halos to examine how the 
background flux from substructure can be used to quan- 
tify the halo shape. We calculate the observed flux as a 
function I and b for spherical, prolate and oblate flattened 
2:1 and 3:1 geometries. In each case, the density profile 
is taken from Eq. |^ and again the observer isplaced in 
either the short (Figure |TT](a) ) or long (Figure fll](b)) axis 
of symmetry. 

These plots show how the distribution of flux on the 
sky can vary significantly depending on the shape of the 
density distribution and on where the observer is situated 
within the halo. 



C. The distribution of point sources 

Individual substructures may be observed and quanti- 
fied if the resolution of the telescope is sufficient. How- 
ever, all of the past and present observations would only 
detect substructure as unresolved point sources. The dis- 
tribution of their fluxes (and spatial distribution on the 
sky) may be used to rule out alternative origins, such as 
extra-galactic sources. In Fig. |l2| we plot the cumulative 
distribution of point sources above a given flux within 
one degree square bins. The two curves consider sub- 
structure with peak circular velocities larger than 10 km 
s _1 and 1 km s _1 . The number density of the brightest 
sources in the sky scales as N oc F~ - 7 . 

Higher resolution simulations are vital to quantify how 
much substructure survives within the galactic halo, how 
it is spatially distributed and to quantify the internal 
structure of surviving substructure. However, Fig. [l^ 
gives an idea of what to expect if an all sky survey is 
carried out that is capable of detecting the brightest sub- 
structure halos. 



IV. CONCLUSIONS 

Numerical simulations that follow the growth of struc- 
ture within a universe dominated by neutralinos (cold 
dark matter) have achieved a resolution that allows their 
global structure and internal structure to be quantified. 
The density profiles, shapes of dark matter halos, abun- 
dance and properties of dark matter substructure, all 
play an important role in determining the absolute sur- 
face brightness of observable products from dark matter 
annihilation. 

We have used the results from the highest resolu- 
tion simulations ever performed of CDM halos to ex- 
amine the expected all-sky distribution of gamma-rays 
from ncutralino annihilation. Substructure can boost 
the expected flux significantly over that originating from 
a smooth dark matter halo. Thus, gamma-ray obser- 
vations, such as EGRET data, may already have the 
potential of constraining a large part of the parameter 
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range of the neutralino cross-sections. The distinguish- 
ing shapes of CDM halos and the unique spatial and flux 
distribution of point sources from substructure within the 
Galactic halo should allow a unique identification of ob- 
servational data with dark matter. 
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FIG. 1. (a) The cir- 

cular velocity curves V c (r) = \J GM(r)/r, and (b) density 
profiles are plotted as a function of the radius for each of the 
halo models considered in the text. 

FIG. 2. The gamma ray flux from neutralino annihilation, 
4>(ip), plotted as a function of the angular distance from the 
galactic center if). The curves show the results using the three 
different density profiles plotted in Fig. ^ The flux at a given 
position is averaged over 47r steradians. 
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FIG. 3. All-sky maps of the gamma ray background con- 
structed using a single high-resolution N-body simulation of 
a cold dark matter halo. The observer has been placed in the 
short (a) and long (b) axis of the simulated halo. 



FIG. 4. The left panel shows a unit oblate ellipsoid and 
the right hand panel shows a unit prolate ellipsoid. The axial 
ratios for both are 2:1. 

FIG. 5. The gamma ray flux, (f>, plotted as a function of 
angle ip, for smooth halos of the same total mass using the 
density profile given in Eq. ^| for spherical, oblate and prolate 
halo geometries. The points are values of the flux measured 
directly from the N-body halo illustrated in Fig. ^[ 

FIG. 6. A sketch showing the geometry of an observer in 
the galaxy viewing substructure in the galactic halo. 

FIG. 7. The gamma ray flux, $av, plotted as a function 
of minimum integration radius R m i n for halo substructure of 
different circular velocities and distances as detailed in the 
text. The shaded region shows the range of background values 
at the Galactic anti-center that can be expected depending on 
the halo shape. The point is the average flux due to all clumps 
with Vpeak > 1 km s . Note that the size of the error bar on 
this point depends on the area of the sky surveyed. 

FIG. 8. The cumulative gamma-ray flux from halo sub- 
structures, Y<4>tot{v > Vp e ak), above a given substructure 
circular velocity v pea k . The ten different curves correspond to 
different Monte-Carlo realizations of a Galactic halo of sub- 
structure halos. The flux is averaged over 47r steradian and 
can be compared with the flux from the smooth halo from 
Fig. | and Fig. § 

FIG. 9. All-sky map of the gamma ray background that 
arises solely from dark matter substructures. The positions 
and circular velocities of sub-halos above a circular velocity 
of lkms -1 are drawn from the N-body simulations but the 
flux from each halo is calculated analytically. The observer is 
located on the short (a) and long (b) axis of symmetry. The 
grey scale corresponds to the log of the flux of annihilation 
products. 

FIG. 10. The average gamma-ray flux per square degree 
from dark matter substructure as measured within the sim- 
ulated CDM halo along a great circle of constant galactic 
latitude (a) and longitude (b). The average has been taken 
over a strip of width 44 degrees, the left hand plot represents 
the view along the short axis while the right hand side is the 
view along the long axis. 



FIG. 11. The effect of halo shape on the gamma-ray flux. 
Halo density profiles are drawn from spherical, oblate or pro- 
late distributions with the indicated axis ratios. The observer 
is placed in the short axis (a) whilst the in (b) the observer 
is in the long axis. 

FIG. 12. The cumulative number of gamma-ray sources 
above a given flux within a window Att = 1° x 1°. The 
two curves are for substructure halos with circular velocities 
larger than lOkms - (dashed line) and lkms -1 (solid line). 
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